Nonlinearity Management in Optics: Experiment, Theory, and Simulation 
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We conduct an experimental investigation of nonlinearity management in optics using femtosecond 
pulses and layered Kerr media consisting of glass and air. By examining the propagation properties 
over several diffraction lengths, we show that wave collapse can be prevented. We corroborate these 
experimental results with numerical simulations of the (2 + l)-dimensional focusing cubic nonlinear 
Schrodinger equation with piecewise constant coefficients and a theoretical analysis of this setting 
using a moment method. 
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Introduction. In the past decade, techniques for man- 
aging dispersion [l| and nonlinearity [3 have attracted 
considerable attention in diverse branches of physics, 
including optics, atomic/condensed-matter physics, and 
other areas. Dispersion management (DM), originally 
proposed for optical fibers based on periodically alter- 
nating the group- velocity dispersion with opposite signs, 
was found to be a robust method for supporting breath- 
ing solitary waves [3- It was subsequently implemented 
in atomic matter waves |^, where the sign of disper- 
sion was controlled via periodic potentials. Nonlinearity 
management (NM) was originally proposed in the con- 
text of layered optical media [a|, but it has garnered 
special attention in the context of ultracold physics Q, 
where it was reformulated as Feshbach Resonance Man- 
agement |6|- In this latter setting, the Feshbach res- 
onance was proposed as a means of periodically mod- 
ifying the scattering length of interatomic interactions 
(and hence the mean-field nonlinearity) to avoid collapse 
in higher-dimensional Bose-Einstein condensates [3] and 
sustain robust breathing coherent structures of atomic 
matter waves [a, Isl . The idea of NM has also inspired 
considerable mathematical research, with a focus on its 
averaged and collapse-preventing properties [ill . 

However, to the best of our knowledge, there has not 
yet been any experiment that addresses the theoretically 
proposed framework of NM. The present Letter aims to 
close this gap by offering an experimental investigation 
of NM in an optical setting. In particular, we study the 
propagation of femtosecond pulses in layered Kerr me- 
dia. In a nonlinear Kerr medium, the index of refraction 
can be written n — no + ri2^, where no is the linear con- 
tribution, n2 is the Kerr coefficient, and / is the beam 
intensity. A beam propagating through such a medium 
with self-focu sing (i.e., 712 > 0) in two-dimensions (2D) 
will collapse [llj if the beam power is above a critical 
threshold, Pc = 7r(0.61)^A^/(8non2), where A is the light 
beam's wavelength. Our medium has layers of glass and 
air, which are both self-focusing but with very different 
Kerr coefficients {n[^^ = 3.2 x IQ-^^cm^/W for glass [l^ 
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FIG. 1: Experimental setup. ND = Neutral density filters, 
NLM = Nonlinear medium, and L1-L3 = lenses. 



and 4^^ = 3.2 x lO-i^cmVW for air [3). The aim 
of our experiment is to showcase what is arguably NM's 
most striking feature — preventing the wave collapse in 
2D that would occur in a homogeneous medium. We il- 
lustrate this effect for powers in the interval (2Pc,6Pc), 
where Pc refers to the critical power in glass, that would 
otherwise be on the verge of collapse before the occur- 
rence of more complex processes, such as the formation 
of a plasma through multiphoton absorption [ij, |l5j. 
The latter offers a defocusing, lossy mechanism that pre- 
vents this catastrophic phenomenon. Optical beams that 
do not diverge due to diffraction (filaments or spatial 
solitons) have been obtained through a variety of phys- 
ical mechanisms |lfil |. Self-guiding and filamentation of 
femtosecond pulses relies on nonlinear losses and neg- 
ative index ch ang es from plasma formation to stabilize 
the beam [ij, [iJ, [iJ, |l^ [l^]. Here we demonstrate a 
new guiding mechanism that does not rely on plasma 
formation and can, in principle, be lossless. The varia- 
tion of the Kerr coefficient in the layered medium pre- 
vents collapse, sustaining an oscillatory variation of the 
beam width for considerable propagation distances (be- 
fore eventual dispersion due to weak losses arising at the 
interfaces between air-glass slides). We compare our ex- 
perimental results to numerical simulations of the non- 
linear Schrodinger (NLS) equation, adapted to the de- 
tailed experiment setup, and find good qualitative agree- 
ment (and quantitative agreement within the appropri- 
ate propagation regime) between the two. An additional 
theoretical understanding of the experimental trends is 
offered by a moment approach Q. 

Experimental Setup. We show a schematic diagram of 




z (mm) 




FIG. 2: (Color online) Dynamics for power P — 5.9Pc- (Top) 
Beam FWHM (in /im) as a function of propagation distance 
(in mm), for air, the layered medium, and glass. The respec- 
tive experimental results are denoted by diamonds, circles, 
and x's (the solid curves are visual guides). The FWHM com- 
puted from Eq. Q is shown by the thin dash-dotted curve, 
and simulations incorporating time-propagation are shown by 
the thick dashed curve. (Bottom right) Spatial dependence 
of the field in Eq. (Q. (Bottom left) The (normalized) ex- 
perimentally measured beam energy for glass and the layered 
medium as a function of the number of glass slides. 



the experimental setup in Fig. ^ A TitaniumiSapphire 
laser amplifier system generates pulses with a duration 
of 160 femtoseconds and an energy of 2 mJ at A = 800 
nm. The beam profile is approximately Gaussian with a 
measured diameter (that is, a full-width half-maximum 
or FWHM) of 5 mm. The beam is attenuated using neu- 
tral density filters (ND) and then focused with a 300 mm 
focal length lens (LI). The nonlinear medium (NLM) is 
composed of 1 mm thick glass slides (Corning MicroSlides 
part number 2947) separated by 1 mm gaps of air. The 
input face of the nonlinear medium is placed 2 mm after 
the beam focus (the beam is diverging when it enters the 
nonlinear medium). At the input face of the nonlinear 
medium, FWHM ss 43/im. After traversing the NLM, 
the beam is attenuated and imaged on a CCD camera 
(Pulnix TM-7EX) using two lenses (L2 and L3) in a 4-F 



configuration, with a magnification of M = 8. We per- 
form multiple experiments in which we vary the number 
of glass slides from one to nine. For each experiment, 
the beam profile is captured both at the output face of 
the NLM and after further propagation through 1 mm of 
air. The 4-F system allows us to image different planes 
along the propagation direction by changing the position 
of the CCD camera with respect to L3. For compari- 
son, we also measure the propagation in glass by placing 
multiple glass slides together (without air gaps). 

Theoretical Setup. The standard model for beam prop- 
agation in optical media, incorporating the dispersive and 
Kerr effects, is the NLS equation [ll|, which we adapt 
to our experimental setting. We rescale space by the 
wavenumber, (■C,??, C) = fc*-^' * (x, y, z) and the electric 
field envelope with u = {n2 /tIq ^^'^E jl3|. (The super- 
script (j) denotes the medium, with j = 1 for glass and 
j = 2 for air.) The physical setting is then described by 

TO,; = -;^V5^u - |upu , < C < [ (glass) , 

1 (1) (2) 

'^C = -o%)^l^-\)\^\'^^ f<C<i(air), (1) 

where glass of scaled length I alternates with air of scaled 
length L — l. Our numerical simulations consider radially 
symmetric solutions of J^l, an approximation supported 
by the experimental data. We compare the simulation 
results directly to the experiments. 

To gain analytical insight, we use the moment 
method [g to reduce the radial NLS equation 
to a set of ordinary differential equations (ODEs) 



for its moments. 



Define /i(C) = J\u\^d^dri, 

4^)(C) = /(|Vw|2-5j>|4)d^dr/ (with gi = 1 and 

92 = ■ni'^ h2^), and /5(C) = \ j \u\'^didri. The quantity 
/i is the beam power and is conserved by the dynamics 
of Eq. (P). The Gaussian profile u = Vbexp{— r^/(2cr^)} 
gives /i = VqIto'^. The remaining quantities are asso- 
ciated with the beam width {I2), momentum (/a), en- 
ergy (/4 ), and nonlinearity (/s). Assuming arg(w) = 
^3''^/(4/2), a good approximation for ^ = in our exper- 
iments, yields a closed set of coupled ODEs for the 7^. 
With the invariants Q^-^'' = 2(4''^ + gjh)^ - ^|/4 and 
Q2 = 2\fTiIr^ l^j, we obtain an Ermakov-Pinney (EP) 
[l9| equation describing the dynamics of the scaled beam 
width y{z) = ^JTM: 
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For the initial Gaussian beam, the invariants are Q\ = 



(2) 



^-a^P\ QX' = iv^''^ /v^'Wi'^ = l-5gi'\ and Q2 = 
TT^a^P^, where a « 1.3556665 and P is the beam's scaled 
(by Pc) power. In the kth segment of the medium, the so- 
-■2 , /-., ^\ 1/2 .., .^ p^ _ r<2 _ 



lution isy= {Ak + BkC + CkC) , with AkBk ~ Cj 
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FIG. 3: (Color online) Comparison of NM for different powers. 
The top left panel shows the experimental results for P = 
2.3Pc (pluses), P = 3.9Pc (stars), P = 4.9Pc (triangles), and 
P = 5.9Pc (circles). We plot the FWHM as a function of 
the propagation distance z. The individual cases of 2.3Pc 
(top right), 3.9Pc (bottom left), and 4.9Pc (bottom right) 
are compared with the PDE diagnostic (see text), shown by 
the dash-dotted curves, and the ODE results (thick dashed 
curves) . 
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FIG. 4: (Color onhne) NM for different layered media. (Top) 
Numerical results for P = 5.9Pc with 1 mm glass-1 mm air 
slide layers (thinnest dash-dotted curve), with 1 mm glass- 
1.5 mm air layers (intermediate thickness dash-dotted curve), 
and with 1 mm glass-2 mm air layers (thickest dash-dotted 
curve). (Bottom) Corresponding solid curves with increas- 
ing thicknesses and symbol sizes represent the experimental 
results for the respective cases. 



1. The initial value j/(C — 0) comes from the experi- 
m.ental setup, and the coefhcients Ak,Bk, and Ck can 
be obtained from continuity conditions at the interfaces 
between the shdes. Results using Eq. |2l will also be 
directly compared with the experiments. 

Results. We applied our experimental, numerical, and 
analytical approaches to a variety of settings. Figure El 
shows the measured beam diameter (FWHM) as a func- 
tion of propagation distance in air, the layered medium. 



and glass. For the layered medium, the beam propagates 
through glass from to 1 mm, through air from 1 mm 
to 2 mm, then through another layer of glass followed by 
air, and so on for nine periods. The pulse power is 5.9 Pc- 
In air, the effect of the nonlinearity is practically negli- 
gible, so the beam diffracts. In glass, the beam focuses 
to a minimum diameter of 7.5 /im after a distance of 2 
mm, defocuses to 14 fim between 3 mm and 5 mm (the 
reason is discussed below), and then starts to diverge. In 
the layered medium, the beam initially focuses through 
the first layer and then its diameter oscillates over the 
first three periods with a mean of about 30 fiia. The 
minimum beam diameter is 20 /xm. The oscillations then 
die down and the beam starts to diverge, although with 
a smaller slope than either propagation in air or in glass. 
Comparing the propagation through the layered non- 
linear medium with the linear propagation (air) , it is clear 
that the beam diameter is sustained over several diffrac- 
tion lengths. In both the glass and the layered medium, 
the beam is self-focused and stabilized. The loss of en- 
ergy is what ultimately causes the beam to diverge. Their 
stabilizing mechanisms are different. The bottom panel 
of Fig.|21shows the total power in the beam as a function 
of the number of glass slides the beam has traversed for 
both the layered medium and glass. In the former, we 
measured a total transmission of 37% for 8 layers (88% 
per glass slide), independent of the power (linear loss). 
This arises mainly from reflections at the air-glass inter- 
faces and perhaps also absorption/scattering. In glass, 
an additional loss appears after 2 mm when the beam 
reaches the minimum diameter. It has been demon- 
strated for similar propagation parameters that the sta- 
bilizing mechanism for propagation in glass is plasma for- 
mation |12| . which occurs through multi-photon absorp- 
tion that creates a negative index change and balances 
the Kerr self-focusing. The multi-photon absorption thus 
accounts for the additional loss measured in glass. The 
resulting larger beam diameter (lower intensity) prevents 
multi-photon absorption from becoming significant in the 
layered medium. The guiding in the layered medium 
is done purely by Kerr self-focusing and diffraction. In 
glass, plasma formation starts at about 2 mm, whereas 
here the beam propagates through the layered medium 
for a much longer distance. We believe the range of sta- 
bilization in the layered medium can be greatly improved 
by reducing the reflection losses due to the refractive in- 
dex mismatch (e.g., by using anti-reflection coating on 
the glass surfaces) . The experimental results are in good 
qualitative and quantitative agreement with the PDE 
simulations of Eq. (Q. In Fig.|21 we show that the beam's 
FWHM accurately follows the propagation in both air 
and glass. In the layered medium, the PDE also follows 
the experiment qualitatively (and even quantitatively at 
first). (Our simulations include the losses at each inter- 
face.) We also used a beam propagation code to perform 
a full 3D simulation that includes dispersion in the tem- 



poral domain. The results are similar to the 2D simula- 
tion, but with an improved quantitative agreement with 
the experiments (Fig. [21 The effect of temporal disper- 
sion is to increase the duration of the pulse, thereby de- 
creasing the power and the strength of the self-focusing. 
In the simulation, the pulse broadens by approximately 
15% after 15 mm of propagation, and the shape changes 
from a Gaussian to a weakly multi-peaked profile. The 
numerical results show that temporal effects do not play 
a critical role in the stabilization of the beam, as shown 
previously for plasma stabilized filaments in glass |l2l | . In 
the rest of the paper, we thus compare the experimental 
results with the 2D simulation. 

We also studied propagation through the layered 
medium as a function of pulse power. The top left panel 
of Fig. 13 shows the evolution of the beam diameter for 
P = 5.9, 4.9, 3.9, 2. 3Pc- We adjust the power by plac- 
ing neutral density filters before the focusing lens LI in 
Fig. n There is a clear trend to stronger self-focusing 
with increasing power. In each case, the beam diverges 
less than for propagation in air. There are some oscilla- 
tions in the beam diameter; they increase with power and 
eventually die down after a few periods as the power de- 
creases. The minimum beam diameter decreases with in- 
creasing power, showing that there is no significant inten- 
sity clamping for these power values. The measured loss 
was the same (within experimental errors) for all power 
levels. The results for 2.3Pc (top right panel), 3.9Pc (bot- 
tom left), and 4.9Pc (bottom right) are also compared 
with PDE and ODE numerical results. The PDF di- 
agnostic (shown by dash-dotted curves) is based on the 
definition of the original beam and correctly follows the 
experiment's qualitative trends (and even its quantitative 
ones for short propagation distances or weaker powers in- 
volving quasi- linear propagation/beam divergence). The 
ODE approximation (thick dashed curves) captures the 
weaker beam divergence for larger P but is less success- 
ful with finer features (such as oscillations in the beam 
width). This may be attributable to the sensitive nature 
of the closure approximations (especially bearing in mind 
the loss properties of the medium). Nevertheless, Eq. Q 
provides a fair, analytically tractable approximation to 
the observations. 

Finally, Fig. 0] shows the effect of changing the peri- 
odicity of the layered structure. The initial power was 
set to 5.9 Pc^ and the air gaps were increased from 1 
mm to 1.5 mm and then 2 mm. The thickness of the 
glass layers was 1 mm in each case. Changing the pe- 
riodicity effectively changes the divergence of the beam. 
As the period is increased, the intensity of the beam is 
lower when it reaches the second and subsequent glass 
layers, resulting in a weaker self-focusing and a "faster" 
divergence. The PDE results for the FWHM once again 
accurately capture the relevant trends qualitatively (and 
also quantitatively at first). 



Conclusions. In sum, we have offered in the present 
work the first experimental realization of nonlinearity 
management (NM) in the context of optical physics us- 
ing femtosecond pulse propagation in layered (glass-air) 
media. We have demonstrated stabilization of the beam 
through NM, which can potentially provide a lossless self- 
guiding mechanism. We have examined the effects of dif- 
ferent beam powers and different layered-media composi- 
tions. We also compared these results with uniform me- 
dia and (partially) accounted for the relevant loss mech- 
anisms. The experimental results are accurately cap- 
tured qualitatively (and, when appropriate, also quanti- 
tatively) by an NLS model and some of their key features 
can also be seen using a far more drastic but analytically 
tractable ODE approximation. Interesting future exper- 
imental directions may involve the reduction of losses or 
the incorporation of slides of defocusing material that 
could lead to a complete stabilization (with stable oscil- 
lations of the beam) |3, Iql ■ 
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